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Abstract. We conjecture on the structure of the quark vacuum from a viewpoint 
somewhat different from, but perhaps supplementary to, standard philosophies. 

Using a rather simple dynamical Hamiltonian model, vacuum excitations carrying 
helicities and ±1 are discussed in connection with the dynamical stability of solutions. 
We speculate on how the rest masses of the light mesons 7r°, ry^and rf could be related 
to these excitations. 
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1. Introduction 

We assume in this paper that the physical vacuum, which has acquired a place of honour 
in contemporary physics because of the great empirical success of the Standard Model, 
is a physical medium that can be probed just like any other physical medium [1]. This 
idea has, of course, a long history behind it. In modern physics, it first emerged as 
a possibihty, suggested by detailed QED studies. It seems, however, that Nature has 
provided us with many alternative ways of empirically testing this idea . We are referring 
in fact to the physics of chiral and flavour symmetries and symmetry breaking in the 
world of quarks and leptons, as assumed and described by the Standard Model [2]. In 
contemporary physics, this is usually tackled by relying on the generally very successful 
methodology of Effective Field Theories [2, 3, 4, 5]. 

In this paper, we concentrate on one of the best hints we have about the nature 
of this presumed medium, through the structure of the lightest self-conjugate mesons. 
This subject has of course been discussed in detail in the last 40 years or so, but we 
believe that the case is not definitely closed as yet. The general question this paper 
studies is: if the vacuum indeed is a mostly unknown physical medium, then how do 
approximate ansatzes about the nature of this presumed medium manifest themselves 
on key properties of the predicted particle spectrum. 

We begin by reformulating this question in terms of Dirac's Quantum Field Theory 
"without dead wood" ([6]; Appendix I). This means that questions about symmetry, 
spontaneous symmetry breaking and chiral anomalies are formulated in a way that 
differ somewhat from the standard one [2, 3, 4, 5]. The main technical points are, first 
we abandon the notion of wave-functionals and vacuum expectation values; this implies 
among other things, in a perturbative setting and at any stage of the development of the 
theory, the automatic non-existence of such irrelevant objects as unconnected vac-to-vac 
diagrams. This non-trivial streamlining of Quantum Field Theory should of course lead 
to the same results as standard formulations, in so far as these agree with experiments. 
However, the interpretation differs. 

We further emphasise that we always work within momentum space. This feature 
we share with pure S-matrix theory, although we have no S-matrix here. 

Specialized principles that we assume are mutually compatible and that we attempt 
to obey are: 

a) essential ideas of the empirically rather successful quark model for mesons [3, 4], 
which are defined only in the so-called confinement phase; 

b) the reliability of the QCD assumptions (especially chiral-fiavour symmetries) for 
the physics at very short distances (i.e. d<^ ^qcd with Aqcd of the order of 0.2 - 0.3 
GeV)[2, 3, 4]; 

c) the presumed moderate reliability of ENJL-fike models for the relevant physics 
of the medium distance quark dynamics (i.e. d > but d< Aq^j^ with of the 
order of IGeV) [7]. We do not actually use this particular model. 

A short preview of this paper would be in order. 
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In Sections 2 and 3 wc present the mathematical machinery used in this and 
subsequent papers. This is based on what Dirac recommended in his later years [6], 
and is eminently suitable to our purposes. It may be rather unfamiliar to many particle 
theoreticians, but the material in App.A may help. 

In Section 4 the master equations previously established are applied to the so-called 
current quarks, which are our fundamental (Weyl) q-numhers (Dirac's terminology, [6]). 
Starting with massless particles in our Model Hamiltonian, we discuss three cases, viz. 
all quarks massless (case A); next , massless (u,d) quarks, massive s-quark (case B); and 
finally, three massive species, but obeying the hierarchy < < <S ms(case C). 

We show that all quarks become massive through pure quantum effects. We study 
the results of self-energy corrections to the quarks energy-momentum relationship, and 
verify the robustness (stability) of the vacuum exhibiting such effects. Particular features 
of this calculation (which to some extent are common to all cases) are: 

(i) automatic output of massive Dirac particles, even with input Weyl fermions; 

(ii) the emergence of Bose-like states (again in Dirac's sense of this word ,Ref 6), 
one of these a massless state (Goldstone,case A), or a nearly massless one (Pseudo- 
Goldstone; cases B and C); 

(iii) the emergence of another massive highly coherent superposition of uu, dd and 

ss. 

The relationship of these vacuum excitations to Tr^,rj^ and rj^' is conjectured upon 
in App.B. 

Wc do not know at the moment how to formulate the problem of confinement in 
our approach (at least as discussed in the literature), let alone find a solution for it. 
We have therefore completely by-passed this issue. All gauge fields of the Standard 
Model (that are assumed to couple to our quarks) provide only some sort of an implicit 
unanalysed background. 

Questions having to do with the electroweak couplings of quarks, although very 
important to our views, shall be discused in a future paper. 

2. The theoretical framework 

The theory used in this paper claims to be yet another example of how to construct a 
"field theory without dead wood" in practice [6] . Some advantages of this more logical 
(and simpler!) approach to field theoretical physical problems (in practice treated by 
standard perturbative methods) were spelled out by Dirac in his later years [6]. 

2.1. The Model Hamiltonian 

We need an effective field theoretic Hamiltonian [2, 3, 4], or a Model Hamiltonian at 
the level of QCD-quarks that must be restricted by the following conditions: 

1) Baryon number and charge conservation ; 

2) CPT-invariance (we further assume CP-invariance) ; 
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3) share with QCD its global colour, flavour and chiral symmetries (we bypass the 
conflnement issue); 

4) have a minimal number of special assumptions and parametrizations. 

The S-matrix at the quark level is not defined, but this effective Hamiltonian should 
ensure that one eventually ends up with an hadronic S-matrix with the usual properties 
[2]. 

There are of course in principle an infinity of "microscopic" quark Hamiltonians 
satisfying these conditions appropriate to the vacuum phenomena discussed in this 
paper, even if we restrict ourselves only to the first two families and to at most 
quadrilinear terms as effective interactions. 

Our model Hamiltonian for QCD-quarks is assumed to have the form: 

H^UB + Hg + Hew (1) 

H^ = H^,, + Vi^^ + Vil^ (2) 

If n represents the quantization box volume, then Ub/^ is a "background energy 
density " as Q ^ oo, supposed to absorb the infinite zero-point energy of the virtual 
Dirac-Weyl quark seas. 

Hoqq, Vqq^ , Vq^ ^nd Hew ^rc picces of our model Hamiltonian trying to capture 
relevant key features of QCD-quarks and their assumed effective interactions as they 
emerge within our model space T>. 

In general we should include in Hq all relevant terms that are at most quadrilinear 
in the preassigned fundamental fermion representations [7]. We conjecture, however , 
that the minimal relevant terms for trying to understand the intricate chiral and fiavour 
dynamics of the lightest quarks can be restricted as follows: 

Hoqqp] ^Y.^^:- ^t^cn - PcpPtp :] (3) 

and 

Vqf [D] ^J:x,[: a+ /3+ + (3,,a,, :] (4) 

C,fJ, 

c,c' /i,ueT> n'yev 

We assume that the input q-numbers etc. are massless colour triplet states 
carrying quantum numbers c and fi = (n.p.X) standing respectively for colour and 
(flavour, 3-momentum, helicity) ; P~ are their antiparticle states, here deflned to mean 
their CP-conjugate states. The symbol :: stands for "normal products" to be presently 
specifled . 

The coefficients are part of the input parameters defining our model space. 
We thus assume that we are dealing with the fiavour representations (with usual 
assignements for u,d,s quarks) 

n = {{l,t3)^{Qi,Si)}® {(0, 0) ® {Qo, So)} (6) 

^ 2 2 
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n — 1 < — > {u, d) 

« = (^,+^)®(|o) d=(i 0) (7) 

n — 2 < — > (s) 

s=(0,0)®(-^,-l) (8) 

The various relevant model matrix elements (5) are specified in Section 5. The 
summation symbols stand for ordinary summations over discrete quantum numbers and 
integrals over 3-momenta. 

It is important to note that interactions among our input quasiparticles are defined 
only within a momentum space T) extending up to about from zero momentum 
(Section 5). 

Colour degrees of freedom are not dynamically involved in this model. We have 
thus nothing to say about confinement as generally understood in the literature. 



3. The approximation scheme 

3.1. Fermions 

The general idea is to introduce a succesion of canonical transformations designed to 
diagonalize as much of this Hamiltonian as possible. The hopefully small and neglected 
residual contributions can be added perturbatively, if necessary. 

Thus, let us begin with the fundamental q-numbers themselves: 

The folowing condition upon the transformation matrix is imposed, ensuring its 
invertibility: 

K\^ + \vi,\^^l (10) 

We furthermore assume (without loss of generality) that these (variational) 
parameters and are constrained to satisfy the conditions 

^ Up, ^ u^^ = COS — v^^vn^v^ = sm— (11) 

hence 

= exp(z^T2) (12) 
where T2 is one of Pauli matrices. 



The quark vacuum 



6 



Our approximation scheme is basically as follows. The q-numbers A^^^ generated 
by the "Dirac-Bogoliubov rotation" (12) should be stationary states (to borrow Dirac's 
terminology [6] , see also Appendix A)) : 

[H,b;]^E,b; (13) 

[H,d;] = E^d; (14) 

We expect that 

H,^Uo + + K^PAb^, b, d+, d) (15) 
H^'^ = T.T. E.KK - d-c-,dt^] (16) 

C jU 

with V^ll being our first "residual" , which is at least quadrilinear. We choose the 
" Dirac-Bogohubov rotation angles" Ofj, so that this residual interaction is as small as 
possible. We should, however, next check whether or not these residuals could destabilize 
our zeroth order solutions. We may then continue with this process, pulling more and 
more out of this V^^} until ordinary perturbation theory eventually may take over, if 
desired. 

In order to reach this goal, we first compute the commutators , linearize the resulting 
equations and look for non-trivial solutions. Using the previously defined effective 
interaction and putting 



Gnxip) 



[ < n'n';X'^\V\nn;Xp> u 



n' \ 
Fnxip) = 



-' n' X' 



where 



(17) 



(18) 



(19) 



fn'X'nx(P,P) =< <^n'p'X''^npX + ^n'fX'f^npX > 

we find Eq.(17) which will be referred to as "the gap equation". Eqs. (18,19) 
represent corrections to the single quark/antiquark energies. 
We find that the condition for non-trivial solutions is: 



and so 



where 



\P\ + Fnxip) - Enx{p) 

en + Gnxip) 



{\P\+F^x{p)}-E^y,{p) 







(20) 

(21) 
(22) 
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Furthermore 

These are the first set of our master equations. Self-consistency requires that non- 
trivial solutions of cqs(20) must be positive definite. 
Solutions are discussed in Section 5. 

3.2. Bosons 

We shall next consider another branch of flavourless boson-like vacuum excitation 
modes. This should give some feeling about the dynamical stability of the above 
vaccum solutions, at least in these particular channels. We continue to work within 
the framework of linearization procedures. 

The Hamiltonian (15) and (16) is thus rewritten in the form: 



H[D] = C/o + E E E / d'pEM{blnp.bcnp^ + dtnfxd-^P^)+ (25) 

c n A 

+ EE E E Jd^^' J d^k<n'^Xs;n' -:^X4V\npXi;n-p\2 



- 2 > 

c.c' ri,n' Ai,A2 A3A4 ' 

We look for stationary solutions of the equation of motion: 

z^6 = [H[V],6] (26) 

A key condition for obtaining such stationary solutions is of course that they should 
obey the conservation laws associated with the symmetries of the effective (model) 
Hamiltonian H[D]. 

We first look for solutions in terms of our primary composites: 

6 = 6(r+,r) 

Tpp' (p) = ^ d-^j_p^jj^fpxa (28) 
In order to ease the notation we conventioned that 

pp' fXaXbip) r etc 
Let us try the ansatz (appropriate momentum integrations understood) 
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r,m r,m r,m 

In this paper, approximations are limited to linearizing the equation of motion. 
This is of course consistent with the basic assumption that our zeroth order solutions 
are truly stable ones, and therefore non-linear corrections are small and expected to 
contribute only to higher orders. So, 

r^(p1 = E / d'^XKMv'nm - E / d'^YKr{p,p')Vrm (30) 
" Jv J. Jv 

with the condition that 

^ / d?mKr? - E / d'^\yKr? = 1 (31) 

^ Jv ^ Jd 

We start the procedure by trying to satisfy the equations 

[^,r+(p-)] = o^(p)r+(p-) (32) 

with the condition that 

^k{p) > (33) 

After some tedious but straightforward calculations we end up with our next set of 
master equations to be used for studying stability questions and conservation laws (in 
short-hand notation): 

Xkpp' = E{~^^/"M'-^M/iVp' +yKt,[i'J^ij.p'pp'} (I) (34) 

Ykp(5' = ^;r~7^^{~^Ktifi'Mp.p'pp' + Xi^fi>J\f^fi>pp,} (II) (35) 

with the definitions: 

PP /AaAfe(p) fifi' n\i\2{k) (36) 

-^pp'M/i' = (37) 
{Ep_ Ep>)5i^p5p'p> + UpiVp < iip\V\ii'p' > u^v^i UpVpi < p'p,'\V\pp, > Up,>Vf^ 

Mpp'p,ji' = 3{m^m^/ < pp'\V\pp' > UpUp>+ (38) 

+VpVp> < p'p\V\ii'fl > v^vl 

J^pp'ij.11' = Hv^Vp' < p'p\V\pp' > UpUp'+ (39) 

+VpVp' < p'p\V\pfi' > UpUpi} 

A comment on effects from couplings (25) in the Hamiltonian [not included in 
Eqs.36 and 37) is in order. Within our linearized theory, it can easily be verified that they 
do not contribute. Their effects, although numerically unimportant within the context 
of our assumptions and approximations, could be treated either by including the terms 
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that were left out and/or by selective summations, as used in standard perturbative 
methods. In the latter case, the Feynman diagram rules must be constructed in such 
a way that one can automatically perform precise corrections for any double countings 
that may arise, due to this early emergence of free point-like bosons in a pure fermionic 
theory. 

It can also be easily verified that this scheme is a conserving approximation, in the 
sense that it satisfies symmetries of the Hamiltonian, and thus conservation laws. Thus, 
for example, let F be the generator of a global gauge symmetry of the Hamiltonian. 
Should the approximate zero-th order vacuum solution not be invariant under F, then 
a zero frequency solution is automatically produced , in agreement with the Goldstone 
theorem (further comments in Section 4). This approximation scheme respects therefore 
the symmetries of the Hamiltonian to leading order of small quantities (i.e. those that 
were left out in the calculations). 

We seek non-trivial solutions for this system of equations. 

4. Numerical Results and Discussions 

We report in this paper the main numerical results based on the theory presented in 
the previous sections. Explicit expressions for the matrix elements (appropriate to our 
model space) are given in this section. 

We assume: 

(i) 

< n'X';n'X';p'\V\nX;nX;p>=< n'n'; X';f/\V\nn; X;p> (40) 

< n'n'; +;0\V\nn; +;p>^< n'n'; +;0\V\nn; -;p >= (41) 

< n'n'; —;;^\V\nn+;p >—< n'n'; —;;^\V\nn; —;p >— — ^ (42) 

and 
(ii) 

< n'+;n'-;:^\V\n+;n-;p>^< n'-;n'+;:^\V\n-;n+;p >^ (43) 

X 

We have thus introduced a minimum of two further (dimensionless) parameters, 
viz. G and gi. Note the explicit - sign on the rhs of definition 42. As we shall see, this 
guarantees positive definite gaps . 

This is supposed to represent effective colour singlet qq interactions at momenta 
below say A-^ fa 1 GeV. This built-in limit thereby acquires a physical significance in 
our model. 
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4.I. Light quark rest masses 

The extensive analysis carried out by Leutwyler et al [5] revealed that the hght quark 

masses obey the condition that their running masses in the MS scheme at scale 
jjL — IGeV must be 

m„ = 5.1 ± 0.9MeV ma = 9.3 ± lAMeV m, = 175 ± 25MeV{U) 
Using eqs. (20-24) and the effective interaction (42) we find that: 

El(x) = el(x) + Al (45) 
en{x) = x- ^vl{x) (46) 
An^xn + gI{An) (47) 
/(An) = E /' x^Un{x)vr^{x)dx (48) 



and 



-2 / \ _ ^ ^n{x) 



= 1 - "'(^) (50) 
with the definitions 

P = xA^ (51) 

An = AnA;, 

en(p) = en{x)A^ 

Enip) = En{x)A^ 

g^l2\G\A^ 

An = en + /(An) 

2) The parameter space must be such that [5] : 

rUsilGeV) = 175 ± 25Mey (52) 

md{lGeV) = 9.3 ± lAMeV (53) 

m^{lGeV) = 5.1 ± 0.9Mel^ (54) 

A sample of our numerical results can be seen in Fig.l through Fig. 8. All fits satisfy 
the additional constraints that |/n| < 10~^ where 

/n = en + /(An)-An (55) 

With these definitions we shall henceforth drop the bars on the symbols in (51). 
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4-2. On the stability of solutions - normal modes 

We shall next examine the dynamical stability of the previous solution. A reasonable 
criterion for stability,or near stability, is that all solutions ftxx of the master equations 
(34 and 35) should be positive or zero, assuming that the gap equation (47) is satisfied. 
The value zero would imply that the chosen vacuum solutions are unstable according to 
this linearized theory. 



4-2.1. Helicity modes These channels are non-trivial solutions of the following 
equations: 

Xafix) = (56) 

u}{x) ^ gjjx) 

_ v){x) ^ _ u){x) ^ 

Yafix) = (57) 

r ( ^^/(^) n ( ^ 

^2E'f{x) + n^ixr"^""' '^^2E'f{x) + n^ixr"^""' 

^2EUx) + n,{xf''^''^ ^2E'f{x) + n,{xf'^''^ 

where 

E'f{x) = 2[Ef{x) - gu}{x)v}{x)] (58) 
and the definitions 

J2 t dx'x'^X^nix, x')ul{x') = A^{x) (59) 

^ /' dx'x'^X,^{x, x')vl{x') = B,{x) (60) 

^ /' dx'x'^Y,n{.X, X')ul{x') = Ca{x) (61) 

J2 t dx'x'^Y,^{x, x')vl{x') = D,{x) (62) 

Putting 

ulix)+vl{x) = l (63) 

^li^) - ^li^) = (^nix) (64) 

K, = ^{A^ + B^-C^-D^) (65) 
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< = |(A,-S, + a-L>,) (66) 
we find that 



This leads to the determinantal equation 

S{uj, x) - 



Sii{uj,x) Si2{u;,x) 

S2l{uj,x) S22{UJ,X)) 



(69) 



where 



5n(.,.)^l-9EP/>V'(^^||l^ (70) 

(71) 

5.(.,.)^2.E^/>^V^p^(|y^;^ (72) 

5.(...) = -l+.i:p/;..V^pgp^ (73) 
where V denotes the principal-value of the integral. 

Again, an essential self-consistency condition for the solutions of these equations 
is that the gap equations must be satisfied. Then it is trivial to show that u = 
( "the Goldstone" ) is a solution, if all parameters that explicitly break the chiral-fiavour 
symmetry of the effective Hamiltonian are set to zero. This is evidently a well-known 
quantum effect. 

Graphical solutions of eq(69) for x —0 can be found from the plots from Fig.l 
through Fig. 8 under various conditions. The zeroes of the function S{uj,0) can be 
interpreted as rest masses of the lowest boson-like vaccum excitations. All but two are 
situated amongst its poles. 

We consider three cases: 

(i) A pure Goldstone mode, where 

Xu^Xd^Xs^O (74) 

In this case, our effective Hamiltonian is chiral - fiavour invariant, but this symmetry 
is broken by the vacuum. A common mass gap in the quark/ ant iquark spectrum thus 
arises out of quantum effects. 

The results in Fig.l and Fig. 2 refer to a toy model, but one that correctly reproduces 
the bare bones of all our calculations, both without self energy corrections and with 
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self-energy corrections. These figures illustrate a simplified case in which all qq input 
energies in Eqs(45) to (50) are of the form 

e(i) = e i = l,2...,Q (75) 

One can then calculate everything analytically. Fig.l illustrates how these mass 
gaps and g are interrelated. We can thus uniquely renormalize g by trading it with the 
in principle "observable" gap. The function S{lj) (Eq. 69) can thus be expressed as 

S{u;) = u^Riu;) (76) 

with R{0) ^ 0. Fig. 2 shows both these functions. Besides the root at zero ("the 
Goldstone") there is another root, roughly 0.7 (arbitrary units). The corresponding 
state vectors can easily be found from the above relevant equations. The heavy state 
is made up of a maximally coherent superposition of all qq basis states. The interplay 
between the single-particle energies and the gap is decisive in determining in which 
frequency range the "heavy" actually shows up. 

The remaining — 1 qq states are orthogonal to this one, and essentially made up 
of a single qq configuration ( i.e. some state i). 

In Fig. 3 the calculations are repeated, but with the input single-particle energy 
spectrum assumed in this paper. 

The quarks energy-momentum relationships within our model space are illustrated 
in Fig. 3 (again, with and without self-energy corrections): no quark species remains 
massless, but the mass degeneracy is of course not lifted. A feature worth noticing is the 
numerical influence of the self-energy corrections on the input single particlc/antiparticle 
spectrum. Again, these arc computed iteratively. About 10 iterations of the gap variable 
are sufficient to reach satisfactory convergence. 

In Fig. 4a and b one can see the relevant branch of vacuum excitations. The quanta 
of these excitations are boson-hke. Note the presence of not only Goldstones (lowest 
rest mass), but also heavy bosons with the largest rest mass, strongly dependent both 
on the q,q spectrum and on mass gaps. Sandwiched in between the poles of 5'(a;), but 
not shown in Figs 4a and b, are other massive excitations with an intermediate degree 
of coherence necessary to make them normalizable and mutually orthogonal. This basic 
structure of our results can be gleaned from Fig.2. 

(ii) A mixed Goldstone- Wigner mode, where 

Xu^Xd^^ ; > (77) 

In this case the quark degeneracy is obviously partially lifted. In Fig. 5 and Fig.6 
one can see the results for a representative g that solves the gap equations. The 
previous massless Goldstone now becomes a "light" pseudo-Goldstone. The move 

towards "constituent masses" for low momenta can still be seen. 

(iii) A mixed hierarquic Goldstone -Wigner mode, where 



a^d > > ; Xs~> x^ 



(78) 
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The relevant results are shown in Figs. 7a,7b and 8. 

Fig 7a shows an overall view of the energies-vs-momenta of the three species of 
quarks (antiquarks). Fig 7b gives a blow-up of this figure for very small momenta. Note 
that the input mass difference — a;^ ~ 0.0042 has become reduced in the output rest 
energy difference E{d, 0) — E{u, 0) ~ 0.00078. This result however does not take into 
account virtual electroweak effects. 

Finally, Fig. 8 shows the corresponding results for bosonic excitation modes. Again, 
the original Goldstone has been metamorphosed into a light pseudo-Goldstone, which 
we now show together with all the heavies. One recalls that our A = basis is in fact a 
linear combination of spin S — and 5" = 1, A = and 1 = and / = 1, /s = states. 
To identify these we need sources which inject these spins and isospins into the vacuum. 
The trends are however unmistakable, even with such relatively unsophisticated effective 
interactions as those adopted here. 

We conjecture that this pseudo-Goldstone would be strongly coupled to the tt" 
whereas 77°' would be strongly coupled to the heavy coherent excitation that was 
previously mentioned (see App.B). The rj^ would correspond to some intermediate 
superposition of the "original" qq background , but where it eventually would land 
strongly depends on secondary residuals in the effective Hamiltonian not considered in 
this paper. 

4-2.2. Helicity ±1 modes The equations determining excitations in these channels are 
formally the same as (56- 58). However, this requires a further parameter (gi), and thus 
we have not pursued the matter in this paper. 

5. Summciry 

That there are self-consistent non-trivial stable solutions of the basic equations (20) 
and (69) suggests that there is indeed support to the conjecture [10] that the quark 
vacuum may well exist in the same universality class as the vacua (ground states) of 
some well-understood superconducting Fermi liquid systems (e.g. '^He in its B-phase, 
[10]). 

If so, we have a well-understood effective field theory that could reaUistically be 
expected to provide support for the standard Chiral Perturbation Theory of pions and 
etas [2, 3, 4, 5] in its attempts to account for features that are not determined by chiral 
symmetry alone. In this connection, we conjecture that Figs 4, 6 and 8 show the rest 
mass spectrum of boson-like particles, where the lowest mass will be strongly coupled 
to 7r°, the next lowest similarly strongly coupled to r^o whereas the heaviest (having a 
maximally coherent superposition of uu, dd and ss) could be the famous t^q. This is a 
robust feature of these calculations. 

The gauge fields of the Standard Model play some sort of a background role not 
analysed in this paper ; neither is the origin of the leading effective interactions defining 
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our model space (see however Ref. [7]), clearly controlling key features of the hghtest 
mesons. 

This conclusion confirms once again Dirac's belief of the fundamental unity of 
Physics [6]. 
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Appendix A. 

We sketch here for the benefit of the reader not familiar with the contents of a "Field 
Theory without dead wood" , the streamhned version of Field Theory proposed by Dirac 
[6], as the most convenient tool for expressing our ideas. 

Given the effective Hamiltonian, our task is to solve the Heisenberg equations of 
motion [6]: 

t^^k{t) = [H,k{t)] (A.i) 

for q- numbers that we presume to be constants of motion K{t). 

We start our approximation scheme by partitioning the Hamiltonian : 

H = H'q + V' (A.2) 

We first diagonalize Hq as best as we can. It may however break some of the 
symmetries of H. Approximation schemes for finding stationary solutions of (A.I) must 
therefore be designed so that they respect the conservation laws associated with these 
symmetries to any desired accuracy. 

Introducing q-numbers in Dirac Picture 

C)^(i) = e'^^*0(t)e-*^°* (A.3) 

T/^(i) = e^^°*y'e-^^°* (A.4) 
we then obtain: 

t^dD(t)^[VD{t),OD(t)] (A.5) 



Setting 



dD{t)^j:d^S\t) (A.6) 

n=0 



we find the conditions: 

^W{t)^0 (A.7) 



dt 
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and the recursion relations 

^^0^S^ it) = [VD{t),dt'\t)] n > 1 (A.8) 

The first condition gives 

d'S>it)^K (A.9) 

where the q-mimbcr K is not cxpUcitly time-dependent. 

The remaining corrections are given by the recursion relation : 

O^D\t) = i-i) f dt,[Vn{h),dt'\t,)] n > 1 (A.IO) 

Before one can interpret the results at any given stage, one must apply the q- 
numbers thus obtained to the standard (Heisenberg) reference ket |0 >t- This reference 
ket is not the usual "vacuum state vector" , as we have no such concept here, but simply 
helps to define normal products at time t only. We adopt the following procedure [6]: 

(i) Compute the commutator and then apply Wick's theorem [8], using of course 
all the previously established relevant conditions, if any. This makes significant the need 
for normal ordering in this q-number theory; 

(ii) Carry out the time integrations; 

(iii) Pick up the coefficient of the i'th term , and interpret this coefficient as the 
prediction/postdiction of this theory for what the amplitude for the i'th configuration 
(state) to emerge at any time t is, i.e.: 

oo oo 

d{t)\o >,= y: oW(i)io >,= y: e{i^ >< ^lo^^'Km >>, (a.ii) 

N=0 N=0 

As emphasized by Dirac, it is not possible within this framework to make 
predictions, or postdictions, of what can be, or could be, observed at any time other 
than t, even if one possesed a complete knowledge of all amplitudes at time t. A critical 
discussion and assesment of this procedure was given by Dirac in the quoted references. 

Appendix B. 

Let us speculate on the structure of the lightest self-conjugate physical mesons at rest. 
We rely on an assertion, common to all versions of the Quark Model, that mesons are 
qq colour singlet composites . So, we define these as the composites 

. JV 

-Y[ d'py.i^I)JxMTi^I)Jx{^) (B.l) 

Putting 



uIJX = M 
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{iq)XiX2 = r ; A = Ai - A2 

we get 
where 

X.^u)jx{p) = E E (^^4 " ^2)1 JA)(i?i - g|70)XM(.,)A,A.(p) (B.3) 

q=—i A1A2 

and mutatis mutandis for the Y amphtudes. 

The amphtudes X^i^nyjxij), p') and Yi,(^ij)jx{p, p')can be determined after finding non- 
trivial solutions of the system of equations (34) and (35). 

Meson-like particles carrying spin J and isospin I at rest would then be describable 
by wave-packets such as 

r+ = / d^p<^M{p)VUp) (B.4) 

where 

Jd'p\^M{p)\'^l (B.5) 

Within these approximations , these composites would appear to be structureless 
point-like spin and spin 1 bona fide free bosons carrying good quantum numbers, in 
all observations involving small momentum transfers.The amphtudes X and Y (B.3) 
will be explicitly needed, should one wish to further test this conjecture, e.g. by 
studying electroweak probings of these point-like mesonic quasiparticles at rest. Some 
residual strong interaction effects, describable within the present scheme, should be 
treated by first returning to our basic theory (Section 2) and going beyond the leading 
approximations, but still within the framework of conserving approximations. In any 
case, this procedure could be a good support to Chiral Perturbation theory as applied 
to pions and etas [3] beyond its leading approximations. 

Specific, but easily handled, field theoretical models such as this one could, among 
other duties, thus help in testing the widespread suspicion that the non-chiral quark 
model description of the lightest mesons, as used in more or less sophisticated so-called 
relativised versions,e.g. [9], is at best too crude to account for these particular mesons 
and their interactions. 
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Figure 1. This shows how the "couphng constant g" is renormahzed in this paper, i.e. 
how it can be uniquely traded with an (in principle) "observable" mass-gap parameter 
A. g is dimensionless, and the gap parameter is in units of 1 GeV. The lower curve a) 
does not include self-energy corrections, but the upper curve b) does. 
The parameters used in this toy model are: e = 0.2 ;f2 = 1000. 



Figure 2. This illustrates how certain relevant vacuum Bose modes arise in the toy 

model : a Goldstonc (playing the role of pions) and a heavy mode (playing the role of 
^0/ 7-j 'pjiig is a highly coherent mode, in contrast to the other modes, wich remain in 
nearly the same state as for 5 = 0, and are sandwiched among the poles of S{uj). The 
parameters used in making Fig2. are: e = 0.2, A = 0.3. 



Figure 3. We plot here the energy - momentum relationship in case A (mass- 
degenerate quarks) . Dimensionsful variables are in units of 1 GcV. No self-energy 
corrections are included in curve E{l,p) , but they are fully included in the curve 
E{2,p). Note the importance (in our case) of the self-energy corrections that tend 
to push the p = masses towards something like "observed constituent " values, 
especially for small momenta. The parameters used in making Fig. 3 and Fig. 4 are: 
A = 0.05 {g = 1.34). 



Figure 4. a) We plot here the function S{u), 0) in case A (massless quarks) . Variables 
having dimensions arc in units of IGcV. We show only the zeros of S{uj) which are 
of interest from the perspective of this paper, i.e. the lightest and the heaviest of the 
excitations. Note the Goldstone at w = 0. 

b) Note the "heaviest", highly coherent excitation at about ui ~ 2GeV. The self- 
consistency conditions are fully respected. Intermediate zeroes have not been shown. 



Figure 5. We plot here the energy - momentum relationship in case B (massless u,d- 
quarks) . Variables with dimensions are in units of 1 GeV. Self-energy corrections are 
fully included. 

The parameters used in making this figure are: Xu = xa = 0; Xg = 0.2 and g = 0.08. 
We find that E{u, d; 0) = 0.008 ; E{s; 0) = 0.205. 
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Figure 6. We plot here the function S{u,0) in case B. Quantities with dimensions 
are in units of 1 GeV. Self-energy corrections are fully included. The same parameters 
as in Fig5. 

Note the disappearance of the Goldstone {S{0, 0) ~ —0.9). Instead, we get a pseudo- 
Goldstone with a rest mass somewhere between w = 0.4 and u = 0.6. 

The "heavy" of special interest to us is the last zero. Its rest mass emerges between 
w = 1.2 andw = 1.4. with our particular choice of parameters. In all cases the details 
of the qq spectrum and the gap parameters are crucial in determining where the masses 
actually emerge. 



Figure 7. a) We plot here the energy - momentum relationship in case C ( 3 massive 
quarks) . The dimmensionfull variables are in units of 1 GeV. Self-energy corrections 
are fully included. The parameters used in this calculation are: 

Xu = 0.0051, xd = 0.0093, Xs = 0.1750 ; g = 0.05 

b) We show in this figure for clarity a blow-up of the low momenta part of a) 



Figure 8. We plot here the function S{u), 0) in case C (massive quarks) . Dimensionsful 

variables are in units of 1 GeV. Only energy corrections on the q,q energies consistent 
with the often mentioned self- consistency conditions are fully included. Same 
parameters as used in Fig. 7a. See text for further comments. 
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